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Shapes and Sizes of Voids in the ACDM Universe: 
Excursion Set Approach 
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ABSTRACT 

We study the global distribution and morphology of dark matter voids in a ACDM 
universe using density fields generated by N-body simulations. Voids are defined as 
isolated regions of the low-density excursion set specified via density thresholds, the 
density thresholds being quantified by the corresponding filling factors, i.e., the frac- 
tion of the total volume in the excursion set. Our work encompasses a systematic 
investigation of the void volume function, the volume fraction in voids, and the fitting 
of voids to corresponding ellipsoids and spheres. We emphasize the relevance of the 
percolation threshold to the void volume statistics of the density field both in the 
high rcdshift, Gaussian random field regime, as well as in the present epoch. By using 
measures such as the Inverse Porosity, we characterize the quality of ellipsoidal fits to 
voids, finding that such fits are a poor representation of the larger voids that dominate 
the volume of the void excursion set. 
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1 INTRODUCTION 

The detailed distribution of mass in the Universe is only 
now being fully revealed by large-scale astronomical sur- 
veys. Consequently, in the last decade or so, a great deal 
of effort has gone into detailed studies of the high den- 
sity regions, in particular the progenitors of galaxies, and 
clusters of galaxies. The density profiles, merging history, 
substructure and many other features have been increas- 
ingly accurately modeled, simulated, and analyzed. This is 
not surprising because these are precisely the places where 
most "action" occurs: the first objects, galaxy formation and 
merging, X-rays, gravitational lensing, etc. 

But even the lower-density, ostensibly "quiet" regions 
of the Universe are of significant interest: recent studies 
indicate a discrepancy between observations and theory in 
low and medium density regions which are naturally associ- 
ated with voids and superclusters respectively. Simulations 
show that dark matter halos in voids are capable of devel- 
oping into void objects observable as (dwarf) galaxies. How- 
ever, the morphology-density correlation so natural in the 
biased galaxy formation scenario seems to be in contrast 
to what is observed. Despite a fairly long history of stud- 
ies [see, e.g., Zel'dovich, Einasto, & Shandarin (1982), Rood 
(1988)] the problems of void phenomenology remain serious 
enough for Peebles (2001) to conclude: "The apparent in- 
consistency between the theory and observations of voids 



is striking enough to be classified as a crisis of the CDM 
model." 

The study of low and medium density regions is de- 
manding from both the observational and modeling perspec- 
tives. From the point of view of observations, large contigu- 
ous regions need to be surveyed in sufficient depth and spec- 
tra, colors, and morphology information for the surveyed 
galaxies obtained. Fortunately, the observational impedi- 
ment is coming to an end: the 2dF redshift survey (Colless et 
al. 2001) has been completed and the data released, and the 
Sloan Digital Sky Survey (SDSS) has begun to release large 
volumes of data (Adelman-McCarthy et al. 2005). As a re- 
sult, valuable information on the luminosity function, colors, 
concentration and other attributes of galaxies in voids and 
void walls have been obtained (Rojas et al. 2004; Hoyle et 
al. 2005). The interpretation of these results, however, relied 
heavily on the purely phenomenological VOIDFINDER al- 
gorithm of Benson et al. (2003) (see also references therein) 
that assumes voids to be nearly spherical and also depends 
on ad hoc parameters such as the criteria for "void" and 
"wall" galaxies first introduced by El- Ad & Piran (1997). 
As a cautionary remark, the relation of the above obser- 
vational results to the physical processes in voids and void 
walls does not seem to be trivial. 

Observational studies of voids and superclusters search- 
ing for various correlations between properties of galaxies 
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and their environment often do not take into account the fact 
that some galaxies may travel distances comparable to the 
sizes of voids (up to roughly half of the radius of a spherical 
void, as estimated in Appendix 1), and are then artificially 
shifted further by the observationally unavoidable mapping 
into redshift space. Therefore the environment where the 
galaxies were originally formed might well be quite different 
from the one where they are ostensibly observed: The physi- 
cal processes that govern the observed properties of galaxies 
in voids and in void walls certainly operate in physical space 
and not in redshift space. 

Additionally, theoretical models of voids remain over- 
simplified; the assumption that voids have close to spherical 
shapes is often the basis of theoretical and observational 
studies (Kauffmann & Fairall f99f; Dubinski et al. 1993; 
Ryden 1995; Amendola, Frieman & Waga 1999; Goldberg & 
Vogeley 2004). Some algorithms formally allow nonspheric- 
ity but in a quite limited form (see e.g., Aikio & Mahonen 
1998). This prejudice is partly based on a result due to 
Icke (1984) who noticed that by changing the sign in the 
solution for the collapse of a uniform ellipsoidal enhance- 
ment (Lin, Mestel & Shu 1965) one can describe the evo- 
lution of an ellipsoidal void. He showed that the aspherici- 
ties of an isolated ellipsoidal underdense perturbation "tend 
to disappear" as the void becomes bigger. Following this, 
many theoretical studies of void evolution have taken the 
spherical model for granted. Icke's conclusion, however, was 
based only on inequalities regarding the relative accelera- 
tions along principal axes of the ellipsoid. He did not inte- 
grate the equations and therefore the evolution of the ra- 
tios of the principal axes remained unknown. Bertschinger 
(1985), who studied the evolution of an isolated uniform 
ellipsoidal void numerically, confirmed the effect found by 
Icke. However, quantitatively his result was not very impres- 
sive: an initially oblate ellipsoid with axis ratio Ci/ai = 0.5 
and density pi — 0.9p(ti) (Si = —0.1) evolved to a more 
spherical ellipsoid with c//a/ = 0.722 when the density in- 
side the ellipsoid dropped to pf = 0.04p(t/) (5f = —0.96) 
after the Einstein-de Sitter universe had expanded to 162 
times its initial size. One can hardly expect that a void can 
evolve for a such long time without being seriously disturbed 
by its neighbors. For instance, as predicted by the adhesion 
approximation, the shape of voids is likely to be affected by 
violent interactions with adjoining superclusters and voids 
resulting in the growth of some voids at the expense of other 
voids (Kofman et al. 1992; Sahni, Sathyaprakash & Shan- 
darin 1994). Another caveat regarding Icke's model is the 
negligence of external shear. Eisenstein & Loeb (1995) have 
emphasized that even in the case of the collapsing ellipsoid 
the shape change is primarily induced by the external shear 
and not by the initial triaxiality of the objects. Due to their 
low internal density, it is clear that the external shear must 
play an even stronger role in the case of voids. In any case, 
the issue of dynamical evolution of voids appears to be un- 
settled and requires further study. 

From the point of view of detailed modeling, voids re- 
quire large simulation volumes, while at the same time, high 
mass resolution is necessary to track the smaller dark matter 
halos that populate the voids. In addition, hydrodynamics, 
feedback, star formation, etc. need to be modeled adequately 
to understand galaxy formation in the void complex. It is 
worth noting here that theory (Kofman et al. 1992; Sahni, 



Sathyaprakash & Shandarin 1994) as well as some N-body 
simulations (Kofman et al. 1992; van de Weygaert & van 
Kampen 1993; Dubinski et al. 1993; Gottlober et al. 2003) 
indicate that both superclusters and voids have nontrivial 
substructure and complicated internal dynamics. As an ex- 
ample, very complex substructures have been detected and 
studied in dark matter superclusters (Shandarin, Sheth & 
Sahni 2004) and voids by the Virgo consortium utilizing N- 
body simulations (Jenkins et al. 1998). 

The problems outlined above require a more systematic 
study of voids - of their morphology and topology, dynamics, 
and internal structure. Assuming that the low dark matter 
density regions are closely related to observational voids, we 
concentrate here on the most basic global properties: their 
sizes and shapes. We study the dark matter mass distribu- 
tion in real space obtained from N-body simulations of the 
ACDM model. Although the mass distribution in real space 
cannot be measured directly, it controls the dynamics, as 
well as other properties of voids. 

It is worth stressing that there is no commonly accepted 
definition of a void. Here we employ the underdense excur- 
sion set approach: voids are defined as individual regions 
where the density contrast is below a certain threshold 5 C . 
The boundaries of the voids are therefore closed surfaces of 
constant density contrast 8 C . The advantages of this defini- 
tion include: (i) simplicity, (ii) a clear uniform definition of 
void boundaries * allowing the study of shapes and sizes, 
and (iii) symmetry with the common definition of super- 
clusters as regions with densities above a specified density 
threshold. We do not preselect any particular threshold but 
rather study the transformation of voids as the threshold 
is raised from low to high values. As in other definitions 
of voids, we also assume a certain smoothing procedure for 
the density field which we take here to be a uniform Gaus- 
sian filter. This allows a continuous transition to a Gaussian 
density field if the scale of the filter window is taken to be 
sufficiently large. 

As further discussed below, one of the consequences of 
this definition is the associated percolation transition: at 
thresholds above a certain value the largest void spans the 
entire volume and quickly becomes the dominant structure 
in the underdense excursion set. At the threshold corre- 
sponding approximately to the percolation transition, the 
individual voids reach their largest sizes and volumes. We 
will quantitatively study the shapes and sizes of underdense 
regions at these thresholds. 

Although we do not attempt to review the entire 
body of work on voids, we would like to emphasize that 
many other definitions of voids have been previously sug- 
gested (Kirshner et al. 1981; Kauffmann & Fairall 1991; 
Sahni, Sathyaprakash & Shandarin 1994; El-Ad & Piran 
1997; Aikio & Mahonen 1998; Plionis & Basilakos 2002; 
Benson et al. 2003; Colberg et al. 2005). All of them arc 
claimed to be in agreement with visual impressions of ob- 
served data and N-body simulations. While "visual impres- 
sion" pertains more to the expertise of psychology rather 
than physics or astronomy, in general these various defini- 
tions agree at least on one basic principle: dynamically, voids 



* Numerous void finding algorithms often define different parts 
of the boundary surface by different conditions. 
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are closely related to the underdense regions. Often the algo- 
rithms finding voids are designed to suppress narrow tunnels 
between two fat voids, see e.g., Kauffmann & Fairall (1991); 
Colberg et al. (2005). This inevitably involves ad hoc param- 
eters describing which tunnel is narrow and which is not. In 
addition, there are attempts to incorporate some concepts of 
theoretical dynamical models into the construction of voids, 
see, e.g., Dubinski et al. (1993); van de Weygaert & van 
Kampen (1993); Benson et al. (2003); Colberg et al. (2005); 
in practice this is essentially the spherical model of void 
evolution. 

In this paper, we develop an unbiased (i.e., without any 
theoretical or aesthetic or other ad hoc assumptions) tech- 
nique that is able to identify voids and then quantify - at 
least crudely - their shapes. The only 'arbitrary' param- 
eter is the scale of smoothing of the dark matter density 
field. The smoothing filter is taken to be a uniform Gaussian 
window [oc exp(— r 2 /2Rj)]. Obviously, this is not an opti- 
mal smoothing method and we are planning to improve this 
part of the technique in future work. A particularly promis- 
ing method of constructing the density field from point-like 
distributions obtained in redshift surveys and N-body sim- 
ulations is based on Delaunay tesselations (Schaap & van 
de Weygaert 2000; van de Weygaert 2002). Density fields 
reconstructed in this manner appear to preserve anisotropic 
features and therefore may have some advantages over con- 
ventional 'Cloud-In-Cell' techniques followed by an isotropic 
smoothing, as performed here. 

As stated earlier, we define voids as isolated regions of 
the low density excursion set specified by a density thresh- 
old: S = Sp/p < S c . This definition is significantly different 
from the many definitions used in the literature in one par- 
ticular sense. It defines the boundary of the void explicitly 
and unambiguously. In contrast, most other studies define 
a void as an underdense region surrounding the minimum 
of the density field (e.g., van de Weygaert & van Kam- 
pen 1993). The boundary may remain unspecified exactly 
(e.g., Gottlober et al. 2003; Sheth & van de Weygaert 2004) 
or be constructed (e.g., Kauffmann & Fairall 1991; Aikio & 
Mahonen 1998). Though in both cases it is implied that the 
boundary approximately corresponds to the surface of con- 
stant density, various exceptions from the rule are usually 
allowed. 

In our systematic study of the void excursion sets at a 
large number of density thresholds, we identify a particular 
threshold when the voids are the largest by volume. This 
threshold is just below the percolation threshold (S < 5 C ) 
where voids merge with each other and form a single per- 
colating void that spans the entire volume of the simula- 
tion (Sheth et al. 2003; Shandarin, Sheth & Sahni 2004). The 
percolating void itself does not have a well-defined volume 
and has a very complicated geometry and topology, how- 
ever, it possesses interesting scaling properties (Bakunin, 
Medvedev & Shandarin 2005) . We show that the individual 
regions of the underdense excursion set do not have spheri- 
cal shapes: the mean ratio of the smallest to largest axis is 
about 0.45 more or less independently of the void volume. 
In addition, using quantitative measures such as the Inverse 
Porosity, we investigate to what extent ellipsoidal and spher- 
ical volumes can be used as good approximations for realistic 
void shapes (as given by our definition). We conclude that 
while such an approach may work for small voids, it is a 



poor approximation for the larger voids which dominate the 
volume of the void excursion set. 

The paper is organized as follows: In §2 we describe the 
N-body simulations used in this study. In §3 we define the 
void volume function, the number of voids of a given vol- 
ume, and the void volume fraction. In §4 we briefly discuss 
voids in a cosmologically relevant Gaussian random field and 
then, in §5, the voids resulting from N-body simulations of 
the ACDM model. In §6 we show how to approximate voids 
by ellipsoids and define the voids' Inverse Porosity. In §7 
we define and quantify the sphericity of voids and discuss 
the statistical distribution of void shapes. In §8 we discuss 
the shapes of the voids in the simulations. Specifically, we 
discuss the mean Sphericity and Inverse Porosity (§8.1) and 
their distribution functions (§8.2) as found from the simula- 
tions. We summarize and conclude in §9. 



2 N-BODY SIMULATIONS 

The dark matter density field at early times is taken to be 
a Gaussian random field, specified completely in real space 
by its two-point statistics and in fc-space, by the primordial 
fluctuation power spectrum P{k), typically assumed to be 
a power-law. At later times, the density field on the length 
scales of interest here is nonlinear and non-Gaussian. Ac- 
curate evolution of the initial density field on these scales 
requires the application of numerical N-body methods. 

Our primary interest here is in the shapes and sizes of 
void regions and not in studying the detailed dynamics and 
galaxy formation within. Since the relevant void sizes are on 
the scales of tens of Mpc, a simulation box of 256 h~ Mpc (h 
is the Hubble parameter in units of 100 km/s/Mpc) provides 
an adequate total volume for the gathering of void statistics. 
A void is defined operationally by considering a smoothed 
density field with a smoothing scale of the order of 1 Mpc; 
it follows that for the studies reported here, force resolution 
demands are modest and easily met. Similarly, as we are not 
interested in the details of the density field interior to the 
void, the mass resolution required is also easily attained. N- 
body simulations with N p = 256 3 particles corresponding to 
an individual particle mass of 1.227-10 11 M are sufficient 
for the task. 

Given the requirements above, a particle-mesh (PM) 
code is completely adequate. The simulations were per- 
formed using the parallel PM code-suite MC 2 (Mesh-based 
Cosmology Code), one of the codes used recently in an 
extensive test of cosmological large-scale structure simu- 
lations (Heitmann et al. 2005). The cosmological param- 
eters were taken to be Q, m — 0.314 (where Q. m includes 
both dark matter and baryons), Q b — 0.044, S7a = 0.686, 
H = 71 km/s/Mpc, ct 8 = 0.84, and n = 0.99. These values 
are in concordance with the WMAP measurement (Spergel 
2003). The matter transfer function was generated using the 
fits from Klypin & Holtzman (1997), the initial redshift was 
z in = 50, and periodic boundary conditions were used. 
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3 VOID VOLUME FUNCTION AND VOLUME 
FRACTION 

The dark matter density field on a grid is obtained from the 
spatial N-body particle distribution by applying the Cloud- 
In-Cell (CIC) method followed by smoothing with a Gaus- 
sian filter. Using this smoothed density field, we define a 
void as an individual region bounded by one or more closed 
surfaces. Thus, a void is identified as a set of grid sites with 
density contrast (6 = (p—p) / p, where p and p axe the density 
and mean density respectively) below a selected threshold 
(S < S c ) and satisfying a neighboring criterion consisting of 
two conditions. First, that the six closest sites are the imme- 
diate neighbors of the selected site and, second, that the im- 
mediate neighbors of every neighbor are also the neighbors 
of the site. The volume of each void is estimated using the 
number of its sites multiplied by the volume associated with 
the grid cell. All voids at a given density threshold make up 
the void excursion set. We quantify the density thresholds 
by the corresponding filling factors (FF), i.e., the fraction 
of the total volume in the excursion set. 

The void volume function n(V) gives the number^ of 
voids having the given volume V 



dN = n(V)dV. 



(1) 



The behavior of the void volume function is shown in the 
left panels of Fig. 1 for a density field at a very early epoch 
z = 150 - essentially a Gaussian field - whereas Figs. 3 and 
4 show the void volume function for evolved nonlinear dark 
matter density fields. These results will be discussed in the 
following sections in detail. 

The volume fraction in voids, v(V) specifies the amount 
of volume in voids of a given volume V, following directly 
from the definition of the void volume function n(V) 

v(V)dV= — n(V)dV, (2) 
Vt 

where Vt is the volume of the simulation box. We stress 
that the very largest void must be excluded from this statis- 
tic because of its peculiar properties around the percolation 
transition, e.g., Shandarin, Sheth & Sahni (2004). We will 
discuss this further below. The void volume fraction is shown 
in the right panels of Figs. 1, 5 and 6 for the Gaussian field 
and N-body simulations. 



4 GAUSSIAN DENSITY FIELD 

We first briefly discuss voids in the case of a Gaussian den- 
sity field. As a cosmologically relevant example of such a field 
we consider the initial particle density field used for the N- 
body simulations. The code generating the initial conditions 
was used to obtain ten realizations of a Gaussian field at the 
stage when ag =< S 2 > x ' 2 = 0.0155 (after smoothing at co- 



moving Rf — 1 h 
field at redshift z - 



_1 Mpc), which corresponds to an initial 
150. 



t It should be noted that in the discussion below, the number of 
voids quoted is for the particular simulation box volume Vb ox = 
256 3 (h^Mpc) 3 and the particular bin sizes: AV/V = 1/V2- 
The statistical significance of the results can be directly assessed 
using the total number of voids in every bin. 





Figure 1. Gaussian field: erg m 0.0155. The panels on the left 
show the void volume function and on the right the void volume 
fraction for different filling factors. The top panels show the both 
functions before the onset of percolation: FF=1%, 3%, 5%, 7% 
and 9% (from bottom to top). The bottom panels show the curves 
after the onset of percolation: FF= 12%, 17% and 25% (from top 
to bottom ). The panels in the middle approximately correspond 
to the percolation transition; the curves correspond to FF C = 10% 
(<5 C = —0.02). The mean and the error bars have been derived 
from ten realizations of the initial conditions. The dashed straight 
line in the left panels show the power law nAV = 8.9x 10 4 V - 0,95 , 
where volume V is measured in (h~ 1 Mpc) 3 and AV/V = 1/vS. 



We first discuss the mean void function from the ten 
simulations. This is displayed in the left panels of Fig. 1. The 
top panel shows the monotonic growth of the void function 
with increase of the void filling factor from FF = 1% (lowest 
curve) to FF = 3%, 5%, 7%, and 9% (top curves) before 
the onset of percolation. The growth stops approximately at 
FF =10% as shown in the middle panel. This corresponds to 
the percolation transition: the merging of voids results in the 
formation of a percolating network in the underdense phase. 
After the onset of percolation the void volume function of 
isolated voids monotonically decreases with increase in the 
filling factor: the curves in the bottom panel correspond to 
FF =12%, 17% and 25% (from top to bottom). 

A robust general scaling relation from percolation stud- 
ies can now be applied to our results. According to this re- 
lation (Stauffer & Aharony 1992), the number of large voids 
scales as 

n(V) xV' T exp(-cV), (3) 

where c tx \FF — FF C | 1/<T ; FF < FF C where FF C is the fill- 
ing factor at the percolation transition, and a and r are real 
parameters. Our results for Gaussian fields are seen to be in 
basic agreement with the percolation model with r « 1.95 
as shown by the dashed straight line in Fig. 1: at percolation 
the void volume function is well-described by a power law 
over three orders of magnitude. Voids with volumes in the 
range of the scaling part of the void volume function would 
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be referred to as "critical clusters" * in percolation theory 
because the power law behavior is characteristic of critical 
phenomena and the properties of these voids are dominated 
by the behavior at FF C . The larger voids in the exponential 
tail of the void volume function [Cf. Eqn. (3)] are rare and 
their properties are not controlled by the system behavior at 
FF C . The scaling properties of the excursion sets in Gaus- 
sian and especially nonlinear fields are interesting topics and 
deserve a more detailed study, however they are beyond the 
scope of the current work. We will return to these issues in 
the future work. 

The volume fraction is shown for the corresponding fill- 
ing factors in the right panels of Fig. 1. The basic behavior 
of the volume fraction is similar to that of the void func- 
tion. The volume fraction grows monotonically as the filling 
factor is increased until the onset of percolation (right top 
panel), after which it monotonically decreases (right bottom 
panel) . The onset of percolation is analogous to a phase tran- 
sition and therefore, due to the expected increase in fluctu- 
ations near a "phase change", it is not surprising that the 
error bars are significantly greater in the vicinity of FF C . 
Although the statistical mean for the percolation transition 
is FF C = 10%, this number will change from realization 
to realization. When analyzing a particular realization, the 
percolation threshold for the particular realization should 
be used rather than the statistical mean value. As a conse- 
quence of this, somewhat different values for the percolation 
threshold than given in this section will be encountered in 
the following analysis. The difference is small, ~ 0.5 — 1%, 
well inside the la range. 



5 VOIDS IN N-BODY SIMULATIONS 

In this section we present a comparative analysis of voids 
in the initial (z = 50) and final (z = 0) density fields 
taken from the N-body simulation. This analysis is targeted 
at establishing which properties of voids are inherited and 
which are acquired in the course of evolution. Figure 2 re- 
lates the filling factor to the density contrast. As we shall 
see later the percolation thresholds at z = incidentally are 
close to 5 = -0.8 (-0.78 and -0.75 for R f = 1 /i _1 Mpc and 
Rf — 2/i _1 Mpc respectively) corresponding to the density 
contrast of the spherically symmetric top-hat void model 
when it reaches the first shell crossing. At these thresholds 
about a half of the volume of the underdense excursion set 
is already in one percolating void. 

It has been demonstrated in previous studies (Shan- 
darin, Sheth & Sahni 2004) that the void volume function is 
dominated by small voids. This is not surprising at the linear 
stage of evolution of the density field, when the field is essen- 
tially Gaussian, but it is also true in the nonlinear regime. 
This is illustrated by Figs. 3 and 4 that show the void vol- 
ume function at different filling factors and smoothing scales. 
In these figures, results are shown for the same underlying 



• In percolation theory, the term "cluster" refers to a group 
of neighboring sites satisfying a criterion such as "occupied" or 
"spin-up", etc. In our case this criterion is set by the density 
threshold <5 C . 
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Figure 2. Void filling factor FF V versus density contrast thresh- 
old & c at z = for two smoothing scales: sml and sm2 corre- 
sponding to Rf = 1 h~ 1 Wlpc and Rf = 2h _1 Mpc respectively. 
The smallest rilling factor is 1% and the largest 35%. 



sml 




Figure 3. Void volume function n(V)AV with the density field 
smoothed at Rf = lh~ x Mpc. The panels on the left are for 
z = and on the right, z = 50. The top panels show the void 
volume functions at small filling factors before percolation occurs 
(FF C =22% , 5 C = -0.78 at z = and FF C = 9.5%, <5 C = -0.14 at 
z = 50) and the bottom panels after percolation. The heavy solid 
lines in the top and bottom panels are the same; they show the 
void volume function approximately at percolation. The dashed 
straight line is the same as in Fig. 1, i.e., the power law nAV = 
8.9 X 10 4 V -0,95 characteristic of the initial Gausssian density 
field. 
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Figure 4. Void volume function n(V)AV with the density field 
smoothed at Rf = 2/i _1 Mpc. The critical parameters are: 
FF C =15% , <5 C = -0.75 at z = and FF C = 9%, <5 C = -0.10 at 
z = 50. The organization of the figure is the same as in Fig. 3. 

density field but with different Gaussian filtering windows, 
Rf = 1 h~ 1 Mpc and Rf — 2 ft _1 Mpc respectively §. 

In Figs. 3-4, the two panels on the left show the number 
of voids at z = while the two panels on the right represent 
the initial density field at z — 50. The top panels show the 
void volume functions at several filling factors correspond- 
ing to relatively low density thresholds when voids do not 
percolate. The bottom panels are complementary to the top 
panels and show the void volume function at the thresh- 
olds after percolation. The heavy solid lines are the same 
in the top and bottom panels and show the void volume 
function approximately at the percolation thresholds. These 
thresholds correspond to a filling factor, FF=9-9.5%, for the 
linear fields regardless of the smoothing scale, and decreases 
from 22% at the smallest smoothing scale (Rf — 1 fo -1 Mpc) 
to 14% at R f = 2/i _1 Mpc for the evolved field. The dashed 
straight line is the power law obtained previously for a Gaus- 
sian field with the same cosmological parameters. 

Figure 3 clearly demonstrates that the evolution results 
in a decrease of the number of small voids and an increase of 
the number of large voids, with the approximate transition 
boundary between the two behaviors at V ~ 300/i -1 Mpc. 
The field smoothed with R f = 2/i _1 Mpc (Fig. 4) yields 
qualitatively similar results, however, there are some quan- 
titative differences. Relative to the Rf = 1 h~ Mpc smooth- 
ing, the number of small voids is lower in both the initial 
and final fields, as well as the number of large voids. The for- 
mer is not unexpected because smoothing with greater scale 
erase smallest voids. However, the latter result seems to be 

§ We also have analyzed the fields smoothed with Rf = 
4/i -1 Mpc but do not show the results here because they do not 
reveal anything significant. 



less obvious since the smoothing scale is relatively small and 
seemed not to influence large voids. Probably the explana- 
tion is in the significant reduction of the percolation filling 
factor in more smoothed density field that makes voids less 
isolated. As a result they merge into the percolating voids 
before they reach large volumes. 

The void volume function at the fully developed non- 
linear stage z = is clearly not a power law, while the 
initial field (z — 50) is mostly in good agreement with the 
scaling discussed in the previous section. The lack of the 
scaling (power law) component in the void volume function 
at the nonlinear stage is a possible indicator of the presence 
of large-scale coherence in the density field. This issue will 
be addressed elsewhere. 

In agreement with the Gaussian case, the results dis- 
played in Figs. 3 and 4 show that the number of voids mono- 
tonically increases with the filling factor and reaches a max- 
imum at the percolation threshold (top panels), then mono- 
tonically decreases with the filling factor (bottom panels). 
We conclude that the largest number of voids - practically 
of every volume - is reached approximately at the perco- 
lation thresholds. Below we show that a small number of 
the largest voids make up most of the volume of the void 
excursion set while the majority of voids have very small 
volumes and constitute only a small fraction of the void ex- 
cursion set. This is true for all filling factors at all smoothing 
scales (Shandarin, Sheth & Sahni 2004). 

Figures 5 and 6 demonstrate - in a manner similar 
to the void volume function (Figs. 3-4) - that the void 
volume fraction grows monotonically with the filling fac- 
tor until it reaches the percolation threshold shown by the 
heavy solid lines (top panels). Thereafter, it monotonically 
decreases with the filling factor (bottom panels). This oc- 
curs because the largest void becomes the percolating void 
and contains most of the volume of the void excursion set. 
The volume fraction of voids reaches a maximum at the per- 
colation thresholds for both linear and nonlinear fields. 

Comparing the void volume fraction (Figs. 5 and 6) 
with the void volume function (Figs. 3 and 4) we find that 
a relatively small number (less than roughly ~ 10%) of the 
largest voids [with approximately V > 10 3 ( /i _1 Mpc) 3 ] make 
up most (roughly ~ 90% or even more) of the void excur- 
sion set (see also Shandarin, Sheth & Sahni 2004). In their 
numbers small voids dominate, but they do not contribute 
significantly to the entire volume of the void excursion set. 

The results of this section are summarized in Fig. 7. The 
top panels show the void volume functions, middle panels, 
the volume fraction and bottom panels, the cumulative vol- 
ume fraction at the percolation thresholds respectively. Solid 
lines correspond to the filtering scale Rf = 1 fo _1 Mpc and 
dashed lines to Rf = 2 ft -1 Mpc. We conclude that the per- 
colation threshold is a natural choice characterized by the 
largest numbers of voids as well as the largest volume frac- 
tion in the void excursion set (compare to Shandarin, Sheth 
& Sahni 2004). As mentioned already, a not completely ob- 
vious result is that smoothing the density field reduces the 
number of large voids along with the expected reduction 
of the number of small voids. It appears that this effect is 
smaller or even absent in Gaussian fields, however, this may 
be a statistical fluctuation due to insufficient statistics. 

The major effect of the nonlinear evolution is sum- 
marized in Fig. 8 where the ratio of the number of voids 
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Figure 5. Void volume fraction u(V)AV with density field 
smoothing of Rf = l/i _1 Mpc, following Fig. 3. 
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Figure 6. Void volume fraction v(V)AV with density smoothing 
scale of Rf = 2h -1 Mpc, following Fig. 4. 



Figure 7. Summary of the information in Figs. 3-6; the volume 
functions and volume fractions are shown at the corresponding 
percolation thresholds. Top: the volume functions at z = (left) 
z = 50 (right). Middle: the corresponding volume fractions. Bot- 
tom: the cumulative void fractions: v dV 
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in the N-body simulations to the number of voids in the 
initial density field at z = 50 is plotted as a function of 
void volume. Voids smaller than 10 3 ( /i _1 Mpc) 3 are strongly 
suppressed while voids with greater volumes become ap- 
proximately three times more abundant than in the parent 
density field. As a result, the voids smaller than roughly 
10 3 ( /i _1 Mpc) 3 comprise less than ~ 10% of the void ex- 
cursion set (the bottom panels of Fig. 7) . Consequently, we 
concentrate mostly on larger voids in our further analysis. 
A somewhat counterintuitive result is that very large voids 



Figure 8. Ratio of number of voids in the N-body simulations to 
the number of voids in the initial density field at z = 50 as a func- 
tion of volume. The voids were identified at the corresponding per- 
colation thresholds in every density field: for N-body simulation 
at FF =22% and 15% at R f = lft- 1 Mpc and R f = 2h" 1 Mpc 
respectively, and for initial state at FF =9.5% and 9% at these 
smoothing scales respectively. 
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[V > 3 x 10 5 ( /i~ 1 Mpc) 3 ] are apparently not more abundant 
at the nonlinear stage than in the initial density field; this 
could be an artifact due to the finite size of the simulation 
box, however, and should be checked using larger volume 
simulations. 

We now discuss how the excursion set definition of voids 
- as regions bound by isodensity surfaces - results in seem- 
ing disagreement with other definitions of voids. Most sig- 
nificantly, we find that the amount of space in isolated voids 
is not large. For a smoothing scale of Rf = 1 h~ 1 Mpc it is 
only about 20% of all space, and for Rf = 2fr~ 1 Mpc it is 
even less: about 15%. In contrast, other studies would claim 
that most of space is taken up by voids. This discrepancy 
is rather easy to explain: We deal with isolated voids nat- 
urally defined by isodensity surfaces, while in other studies 
voids are built using procedures that follow the density con- 
tours only approximately, thereby allowing closure of the 
boundaries of the constructed voids by non-isodensity sur- 
faces. With our definition of voids one can get as much of 
void space as one wants by raising the density threshold, 
but then the natural isodensity surfaces would inevitably 
build a single void percolating throughout the entire volume 
plus a few small isolated voids. Cutting the percolating void 
into pieces would finish the job of constructing large "iso- 
lated" voids, however this brings a certain arbitrariness, and 
therefore bias, into the analysis. Other algorithms actually 
suggest various methods of doing this, but without explicit 
discussion. Our procedure does not cut the largest void - nei- 
ther before nor after it percolates - and therefore the voids 
found by it can make up only a relatively small fraction 
of space. At even greater density thresholds the percolat- 
ing void is absolutely dominant (Shandarin, Sheth & Sahni 
2004). Reducing the smoothing scale increases the number 
of voids and the amount of space occupied by isolated voids 
(Fig. 7). However, if a uniform filter is used for smoothing 
then the intrinsic discreetness of N-body simulations takes 
over and prevents further reduction of the smoothing scale. 
An adaptive smoothing filter would help to ameliorate this 
problem. 



6 FITTING TO ELLIPSOIDS 

As discussed in the Introduction, the actual shapes of voids 
are often taken to be approximately spherical or ellipsoidal. 
In order to investigate possible systematic difficulties with 
this kind of fitting, we identify all voids at a selected density 
threshold (roughly at the percolation threshold) and then fit 
each void by an ellipsoid having the same inertia tensor as 
the void itself: 

J xx = ^2 m i(Vi + Z i )' Jxy = miXi V^ ( 4 ) 

i i 

plus cyclic permutations, 

where rrii and Xi,yi, Zi are the mass and coordinates of the 
site with respect to the void center of mass and the sum is 
taken over all sites belonging to the void. The semiaxes of 
the fitting ellipse can be found from the principal moments 
of inertia Ji, J2, J3, viz., 

a = —rj(J2 + Js — Ji) plus cyclic permutations, (5) 



where M = rrii is the mass of the void. Our focus here 
is on the geometrical properties of voids and therefore we 
compute the inertia tensor of an empty void by setting all 
the rrii to a constant. We reserve the study of another inter- 
esting option, where the rrn are actual dark matter masses 
in the grid cells, for future work. Figure 9 illustrates the 
fitting procedure outlined above for two example voids: in 
the figure, the voids are shown on the left with the cor- 
responding ellipsoids on the right. The axes of the fitting 
ellipsoid, a > b > c, can be used to evaluate the volume of 
the ellipsoid Ve = (4n/3)abc as well as the ratios b/a and 
c/b which approximately characterize the shape of the void. 
The Planarity P = (1 — b/a) /(l + b/a) and Filamentarity 
F — (1 — c/b) /{I + c/b) (Sahni, Sathyaprakash & Shandarin 
1998; Shandarin, Sheth & Sahni 2004) are equivalent mea- 
sures of void shape. 

If a void is ideally ellipsoidal then its volume must equal 
that of the fitting ellipsoid. Actual voids are not ellipsoidal in 
shape, however (Fig. 9), thus their volumes Vv are always 
smaller then the volumes of the fitting ellipsoids Ve- The 
ratio P — Ve/Vv - dubbed the Porosity - quantifies the 
quality of the ellipsoidal fit. We use the Inverse Porosity 
IP — Vv /Ve (0 < IP < 1) as a simple estimate of the 
goodness of the ellipsoidal fit. If IP is close to unity, the fit 
may be regarded as good, otherwise the smaller the IP, the 
worse the fit. 



7 SPHERICITY 

Deforming a sphere into an ellipsoid with three axes, a > 
b > c, one obtains a two-parameter family of shapes that 
can be fully characterized by two ratios: b/a and c/b. Never- 
theless, we would like to approximately describe the depar- 
ture from a sphere with only one parameter. Ideally the pa- 
rameter should have a simple geometric interpretation and 
only a small variation when the remaining degree of free- 
dom is varied. We suggest the use of the ratio of the small- 
est to largest axis, a — c/a, for this purpose. Using this, 
one can distinguish three archetypal shapes among the en- 
tire two-parameter family characterized by a: pancake-like 
with b/a = 1 and c/b — a, filament-like with b/a = a and 
c/b — 1, and ribbon-like with b/a — c/b — \fa.. All three 
have very distinct limits at a — > 0. Figure 10 provides a vi- 
sual illustration of the variation of shapes from a pancake 
through a ribbon to a filament for three values of a. 



8 SHAPES OF VOIDS IN N-BODY 
SIMULATIONS 

In this section we discuss the shapes of voids in the simu- 
lated ACDM model at the present epoch (2 = 0). As earlier, 
we analyze the density distribution in real space at three 
different smoothing scales: l/i _1 ,2/i _1 , and 4ftr 1 Mpc. 

8.1 Mean values 

The mean Sphericity as a function of the volume is shown 
in Fig. 11. Every line shows the ratio at a different filling 
factor and for all three smoothing scales. The variance is 
shown by vertical bars for one of the lines, the others being 



© 0000 RAS, MNRAS 000, 000-000 



Shapes and Sizes of Voids 9 




Figure 9. Comparison of an ellipsoid fit to high and low porosity voids. The low porosity small void (bottom left) has volume of 
V = 7500 cells, Porosity P = V E /V V = 1.92, and Inverse Porosity IP = V V /V E = 0.52; the fitting ellipsoid is at the bottom right. The 
high porosity large void (top left) has volume of V = 24000 cells, P = 4.72, and IP = 0.21. 



very similar. It is clear that there is no significant statisti- 
cal difference between the mean shapes of voids at different 
filling factors or smoothing scale: on average b/a w 0.65 and 
c/a w 0.45. Thus, if a region of the underdense excursion 
set is approximated by an ellipsoid then the ellipsoid should 
be triaxial with the shortest axes at least twice shorter than 
the longest one. A substantial number of ellipsoids has a 
shortest axes as small as a third of the longest one. Figure 
10 provides the visual illustration of such ellipsoids. 

In Fig. 12, we illustrate the dependence of the mean 
Inverse Porosity IP — Vv /Ve as a function of the volume 
of the void for different filling factors and smoothing scales. 
Following the behavior of the Sphericity, the Porosity does 
not show a strong dependence on the filling factor, however 
there are two other characteristics that are very different 
from the mean Sphericity. The first is a systematic growth 
of the mean Inverse Porosity with the smoothing length, 
albeit with large variance. This is an obvious consequence 
of smoothing; the voids in the smoothed density fields being 
more ellipsoidal. The other is a clear dependence on the 
volume of the voids: the larger the volume, the smaller the 
Inverse Porosity and therefore the worse the ellipsoidal fit. 
Combining the information in Figs. 12 and 7 we conclude 
that large voids V > 3 x 10 3 /i _1 Mpc 3 - making most of 
the void excursion set - cannot generically be adequately 
approximated by ellipsoids. 



8.2 Distribution Functions 

Statistical information regarding the quality of elliptical and 
spherical void fits can be gleaned from the distribution func- 
tions for the Inverse Porosity and the Sphericity. We begin 
by first discussing the Inverse Porosity distribution func- 
tion. Figures 13 and 14 show the IP distribution for density 
fields smoothed on three scales: Rf = 1 h~ Mpc (solid line), 
Rf = 2ft~ 1 Mpc (short-dashed line), and Rf — 4/i _1 Mpc 
(long-dashed line). Figure 13 includes all the voids with 
volume greater than 10( /i _1 Mpc) 3 and Fig. 14 those with 
volume greater than 1000( /i _1 Mpc) 3 . The first cut is im- 
posed because the shape measurement of small voids V < 
10( /i _1 Mpc) 3 is not reliable. The second cut selects only 
large voids V > 1000( /i _1 Mpc) 3 comprising approximately 
90% of the void excursion set (Cf. Fig. 7). If voids with 
V < 1000(/i _1 Mpc) 3 are included, the IP distribution peaks 
relatively close to unity, implying that elliptical fits are quite 
good for the majority of the voids. As we demonstrated ear- 
lier, however, small voids dominate only in number and not 
in volume. If small voids are excluded (Fig. 14) then the IP 
distribution peaks at Vv /Ve ~ 0.45 for Rf — 1 ft _1 Mpc and 
Rf — 2/!T 1 Mpc. In this case the ellipsoidal fit is quite poor 
for most of the larger voids that make up roughly 90% of the 
void excursion set. Therefore, in the following analysis we 
will exclude voids with volumes less than 1000( /i _1 Mpc) 3 . 

We note that the IP distribution in Fig. 14 displays 
a trend of the distribution peak shift to the right with 
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Figure 10. Three principal projections of an ellipse arc shown 
in every panel. Each column shows one of three archetypes of the 
deformed sphere parameterized by one number ct = c/a(0<a< 
1): pancake-like (b/a = 1, c/b = a), ribbon-like (b/a = ^/a, c/b = 
\/a), and filament-like (b/a = a, c/b = 1). 



Figure 11. Axis ratios as a function of volume. The average of 
the axis ratios is independent of the filling factor and smoothing 
length. Except for very small voids, the average is constant with 
volume. 



increased smoothing scale, with an especially pronounced 
change for Rf = 4/iT 1 Mpc. This is expected as smoothing 
will improve the ellipsoidal fit for the smallest voids consid- 
ered; for Rf — 4/i~ 1 Mpc, a sphere with the radius of the 
smoothing length has a volume of 270( /i _1 Mpc) 3 which is 
more than one quarter of the volume of the smallest voids 
that dominate the void number distribution. 

Turning to the Sphericity, Fig. 15 shows the Sphericity 
distribution function for large voids Vv > 1000( ft _1 Mpc) 3 . 
In this case, the results from fields smoothed on different 
scales are not very different. The mean Sphericity of the fit- 
ting ellipsoids is about 0.45 (in agreement with Fig. 11). The 
middle panels of Fig. 10 show three archetypal ellipsoids cor- 
responding to this value. The top and bottom panels show 
the range of Sphericity for most of the voids. 



9 SUMMARY AND DISCUSSION 

In this paper, we studied the distribution of underdense re- 
gions (voids) in N-body simulations for the ACDM model, 
voids being defined as individual regions of the underdense 
excursion set: S < S c . In practice, we additionally applied 
the neighbor of neighbor criterion for grid sites satisfying 
the above threshold condition. By defining the void volume 
function, we then showed that the largest number of voids 
of any size is reached just before the onset of percolation. In 
other words, as the density threshold is increased more and 
more voids are found and the total volume fraction in voids 
increases; this continues up to a critical point beyond which 
the number of voids and the volume fraction in all but the 
percolating void begin to decrease. 

The percolation transition motivates a particular den- 
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Figure 12. Inverse Porosity (IP = V/Vt) as a function of void 
volume. The larger the volume the larger is the Porosity on av- 
erage. There is little dependence on the filling factor, while the 
smoothing length systematically suppresses the Porosity. 



sity threshold and corresponding filling factor to focus fur- 
ther analysis, i.e., the percolation threshold itself. After per- 
colation, though the volume of the underdense excursion 
set continues monotonically to increase with the growth of 
the density threshold, the number and volume of individual 
voids decrease because merging of individual voids into the 
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Figure 13. Normalized distribution functions of the Inverse 
Porosity (Vy- /Ve) of voids with volumes Vy > 10( /i _1 Mpc) 3 
at three smoothing scales. The distribution functions arc shown 
for the corresponding percolation thresholds. 



Figure 15. Sphericity distribution function of voids with vol- 
umes Vv > 1000( /i _1 Mpc) 3 at three smoothing scales. The dis- 
tribution functions are shown for the corresponding percolation 
thresholds. 




0.2 0.4 0.6 0.8 

Inverse Porosity 



Figure 14. Normalized distribution function of the Inverse 
Porosity (V v /V E ) of voids with volumes V v > 1000( /i^Mpc) 3 
at three smoothing scales following Fig. 13. 

single percolating void dominates. The percolating void con- 
tinues to grow monotonically with the threshold, mostly by 
absorbing individual voids. Choosing a particular threshold 
above the percolation point seems to be arbitrary in most 
cases and may significantly affect the results one obtains. 

One could cut the percolating void found at an arbitrary 
threshold into pieces. But without any clear physical motiva- 



tion this would introduce an additional arbitrariness which 
we wished to avoid. It should also be noted that special pre- 
caution should be exercised in one's notion of a void, as it 
is always possible to design an algorithm producing more 
or less spherical voids if so desired. The guideline provided 
by Icke (1984) suggesting that voids become more spherical 
in the course of dynamical evolution should be reconsidered 
in a more realistic environment rather than in the case of 
a single void placed in a homogeneous universe. The exter- 
nal shear is likely to play an important role (Eisenstein & 
Loeb 1995) and must be included in the model. Even in the 
case of a single void the effect may be quantitatively insuf- 
ficient (Bertschinger 1985). We will address these issues in 
future work. 

As a cautionary note we stress that the reported re- 
sults may depend on the smoothing procedure adopted here 
and the values of the filtering scale. This is an issue with 
essentially all proposed definitions of voids and needs to be 
further investigated. 

The density field at an early enough epoch is essentially 
Gaussian random on scales below the smoothing scales em- 
ployed here. We considered ten realizations of a density field 
at z = 150 and showed that the behavior of the void volume 
function closely followed a mathematical form motivated by 
percolation theory; at later redshifts, the field is no longer 
Gaussian and this form does not hold. 

At the present epoch, we found that less than 10% or 
so of the largest voids, with volumes greater than about 
10 3 ( /i _1 Mpc) 3 , hold more than 90% of the void excursion 
set while small voids dominate in their numbers. The dy- 
namic evolution of voids results in a significant decrease of 
the number of small voids [V < 10 3 ( /i _1 Mpc) 3 ] and an in- 
crease in the number of large voids. 

In order to study effective void shapes, we fit the un- 
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derdense regions to ellipsoids with the same inertia tensor 
as the voids. We used this to define the Sphericity and (In- 
verse) Porosity of the voids and showed that most large voids 
are not spherical and have irregular geometries which lead to 
high porosity (low Inverse Porosity). The prevailing assump- 
tion of spherical voids, used extensively in the literature, is 
not a good approximation if the actual isodensity surfaces 
are employed as the void boundaries. We found that the ra- 
tios of the semi-major axes have mean values of 0.65 and 
0.45 and the Sphericity has a mean value of c/a w 0.45, 
(where c and a are the smallest and the largest semiaxes 
of the fitting ellipsoids), which is clearly not spherical. In 
addition, the Porosity measure showed that large voids are 
quite porous whereas small voids are more regular, with even 
ellipsoidal fitting becoming problematic for the larger voids. 

Two questions are of potential relevance here. First, the 
representing ellipsoid has the same inertia tensor as a homo- 
geneous void while actual voids are inhomogeneous. There- 
fore one could ask: would the true density profile change the 
fitting ellipsoid to be closer to a sphere? The obvious answer 
is no; it is natural to expect the central part of a void to be 
more spherical than the boundary region but this would have 
a smaller effect on the inertia tensor since the density at the 
center of the void is lower than near the boundary. 

The second question relates to the shapes of voids before 
the onset of percolation. Is it possible that the percolating 
void has the shape of a set of spheres connected by narrow 
tunnels - one of the major concerns of Kauffmann & Fairall 
(1991) and El-Ad & Piran (1997)? This appears unlikely be- 
cause the percolating void is made by the merger of individ- 
ual voids as the density threshold is increased. The merger 
process happens quite quickly near the threshold (i.e., with a 
small change of the threshold) therefore the individual voids 
making up the percolating void cannot change their shapes 
appreciably. Our results show that the largest of them are 
neither spherical nor even elliptical. Additionally, study of 
individual voids before and after percolation do not show 
much difference (Shandarin, Sheth & Sahni 2004). 
Acknowledgments 

S.H. and K.H. acknowledge support from the Depart- 
ment of Energy via the LDRD program of Los Alamos Na- 
tional Laboratory. The calculations described herein were 
performed primarily using the computational resources of 
Los Alamos National Laboratory. A special acknowledgment 
is due to supercomputing time awarded to us under the 
LANL Institutional Computing Initiative. This research is 
supported by the Department of Energy, under contract W- 
7405-ENG-36. H.A.F. was supported in part by the GRF at 
the University of Kansas. S.S. was partly supported by NSF- 
RP087 grant. S.S., K.H., and S.H. acknowledge support from 
the Aspen Center for Physics where the manuscript was fi- 
nalized. We acknowledge useful discussions with Ravi Sheth 
and especially Rein van de Weygaert whose detailed com- 
ments were invaluable. 



REFERENCES 

Adelman-McCarthy, J.K., et ai. 2005, astro-ph/0507711 
Aikio, J., Mahonen, P. 1998, ApJ, 497, 534 

Amcndola, L., Fricman, J. A., Waga, I. 1999, MNRAS, 309, 465 
Benson, A.J., Hoyie, F., Torres, R., Vogeley, M. 2003, MNRAS, 
340, 160 



Bakunin, O., Medvedev, M.M. & Shandarin, S.F. 2005, to be 

published 
Bertschinger, E. 1985, ApJS, 58, 1 

Bond, J.R., Kofman, L. & Pogosyan, D. 1996, Nature, 380, 603 
Colberg, J.M., Sheth, R.K., Diaferio, A., Gao, L., Yoshida, N. 

2005, MNRAS, 360, 216 
Colless M.M. et al. 2001, MNRAS, 328, 1039 

Dubinski, J., da Costa, L.N., Goldwirth, D.S., Lecar, M., Piran, 
T. 1993, ApJ, 410, 458 

Eisenstein, D.J., Loeb, A. 1995, ApJ, 439, 520 

El-Ad, H., Piran, T. 1997, ApJ, 491, 421 

Freedman, W.L., et al. 2001, ApJ, 553, 47 

Goldberg, D.M., Vogeley, M.S., 2004, ApJ, 605, 1 

Gottlober, S., Lokas, E.L., Klypin, A., Hoffman, Y. 2003, MN- 
RAS, 344, 715 

Grogin, N.A., Geller, M.J. 2000, AJ, 119, 32 

Heitmann, K., Rickcr, P.M., Warren, M., Habib, S. 2005, ApJS, 
160, 28 

Hoyle, F., Vogeley, M.S. 2004, ApJ, 607, 751 

Hoyle, F., Rojas, R.R., Vogeley, M.S., Brinkmann, J., 2005, ApJ, 

620, 618 
Icke, V. 1984, MNRAS, 206, lp 

Jenkins, A.R., Frenk, C. S., Pearce, F. R., Thomas, P. A., Colberg, 
J. M., White, S. D. M., Couchman, H. M. P., Peacock, J. A., 
Efstathiou, G., Nelson, A. H. 1998, ApJ, 499, 20 

Kauffmann, G., Fairall, A.P. 1991, MNRAS, 248, 313 

Kirshner, R.P., Oemler, A., Schechter, P.L., Shcctman, S.A., 1981, 
ApJ, 248, L57 

Klypin, A.A., Holtzman, J. 1997, astro-ph/9712217 

Kofman, LA., Pogosyan, D., Shandarin, S.F., Mellot, A.L. 1992, 
ApJ, 393, 437 

Lin, C.C., Mestel, L., Shu, F. 1965, ApJ, 142, 1431 

Peebles, P.J.E. 2001, ApJ, 557, 495-504. 

Peebles, P.J.E. 1980, The Largc-Scale Structure of the Universe, 

Princeton University Press. 
Plionis, M., Basilakos, S. 2002, MNRAS, 330, 399 
Rojas, R.R., Vogeley, M., F. Hoyle, Brinkmann, J. 2004, ApJ, 

617, 50 

Rood, H.J. 1988, Ann. Rev. Astron. Astrophys. 26, 245 
Rydcn, B.S., 1995, ApJ, 452, 25 

S.F., MNRAS, 282, 641 
Sahni, V., Sathyaprakash, B.S. Shandarin, S.F. 1994, ApJ, 431, 

20 

Sahni, V., Sathyaprakash, B.S. Shandarin, S.F. 1998, ApJ, 495, 
L5 

Schaap, W.E. and van dc Weygaret, R. 2000, A&A, 363, L29 
Schmidt, J.D., Ryden, B.S., Melott, A.L. 2001, ApJ, 546, 606 
Shandarin, S. F., & Zcl'dovich, Ya. B. 1989, Rev. Mod. Phys., 61, 
185 

Shandarin, S.F., Sheth, R. & Sahni, V. 2004, MNRAS, 353, 162 
Sheth, J., Sahni, V., Shandarin, S.F., Sathyaprakash, B. 2003, 

MNRAS, 343, 22 
Sheth, R.K., van de Weygaert, R. 2004, MNRAS, 350, 517 
Slezak, E., de Lapparent, V. & Bijaoui, A., 1993, ApJ, 409, 517 
Spergel, D.N, et al. 2003, ApJS, 148. 

Stauffcr, D., Aharony, A. 1992, Introduction to Percolation The- 
ory, Taylor & Francis, London, Washingon, DC 

van de Weygaert, R. 2002, "Froth across the Uninverse, Dynamics 
and Stochastic Geometry of the Cosmic Foam" , arXiv: astro- 
ph/0206427 

van de Weygaert, R., van Kampen, E., 1993, MNRAS, 263, 481 
Zel'dovich, Ya.B., Einasto, J., Shandarin, S.F. 1982, Nature, 300, 
407 



© 0000 RAS, MNRAS 000, 000-000 



Shapes and Sizes of Voids 13 



APPENDIX A: ESTIMATE OF MATTER FLOW 
IN VOID FORMATION 

We suggest the following simple estimate of the comoving 
distance traveled by a fluid element in the process of the 
formation of a spherical void of radius R. We assume: (i) 
the initial density is the mean density, (ii) the final density 
is uniform in the void and is a fraction / < 1 of the mean, 
(iii) the motion is purely radial, and (iv) all the mass is 
accumulated on an infinitesimally thin shell of radius R. 

The model is purely kinematic and makes no assump- 
tions regarding the dynamics. It corresponds to the following 
mapping: 



q + Aq 

R 



if 
if 



q + Aq < R, 
q + Aq > R, 



(q < q c ), 
(q > qc), 



(Al) 



where q and r are the initial and final positions of a fluid 
element. The amplitude A is determined by the value of /: 
A = J -1 / 3 — 1. The critical value q c where the mapping 
changes its form is q c = R/(l + A) = f 1 ^ 3 R. 

The mean distance traveled by a fluid element is 



AtyR 3 
f 



q c ±1 

J (Aq)A%q 2 dq + J (R-q)4nq 2 dq 



1 



R, 



(A2) 



(A3) 



which obviously understates the actual mean distance due 
to assumption (iv). The maximum distance traveled by fluid 
elements is d max = (1 - f 1/3 )R. If / = 0.1 then d « 0.22 
and dmax ~ 0.54. 

Denoting the fraction of the distance by £ = d/R one 
can obtain the pdf of the distances traveled by the fluid 
elements 



dF = 3 



l-2£ + 

o < C < i - f 1/s 



i + 



/ 



(1-/1/3) 



(A4) 
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